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1.0  INTRODUCTION 


The  parallel  plate  waveguide  consisting  of  two  perfectly  conducting  plates  enclosing 
a  homogeneous  isotropic  medium  is  very  well  known.  When  it  is  assumed  that  the  z 
direction  is  the  direction  of  propagation  and  the  plates  are  located  along  the  x  direction 
(atx  =  0  and  x  =  a)  (see  Figure  1),  the  mode  solutions  along  the  x  direction  are  sinusoidal, 


either  sin 


f  nrcx  ^ 


or  cos 


(  nnx^ 


,  where  a  is  the  plate  separation  and  n  can  be  any  integer. 


V  a  )  \  a  ) 

These  sinusoidal  functions  are  periodic,  and  higher  order  modes  have  nodes  that  can 
segment  the  plate  separation  distance  into  equidistant  parts  (Reference  1).  Since  n  can  be 
any  integer,  the  number  of  such  modes  at  a  given  frequency  can  be  denumerably  infinite. 


In  this  report,  a  parallel  plate  waveguide  containing  an  inhomogeneous  dielectric  (or 
pennittivity)  function  of  specific  form  is  considered.  It  is  found  that  by  assuming  an 
exponential  dielectric  function  of  position,  it  is  possible  to  remove  one  of  the  waveguide 
conducting  plates  and  still  keep  the  fields  closely  confined  to  the  remaining  plate  and 
held  within  the  inhomogeneous  region. 


2.0  INHOMOGENEOUS  PERMITTIVITY 
IN  MAXWELL’S  EQUATIONS 

Assuming  that  the  propagation  direction  is  along  the  z-axis  and  that  the  z-direction 
dependence  and  the  time  dependence  are  given  by  ;  We  start  with  Maxwell’s 

equations  in  the  form 

V  x  E  =  icoB\  VD8  =  0 

and  (1) 

VxH  =  -icoD;  VDD  =  0 

Assuming  that  we  have  a  permittivity  function  depending  on  a  single  transverse 
coordinate  and  that  the  permeability  is  constant,  let 

B  =  ju0/uH  (2a) 


and 


D  =  s0s[x)E 


(2b) 
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Furthermore,  assuming  an  infinite  flat  slab  type  of  geometry  (see  Figure  1)  with 

d 

—  =  0 ,  Maxwell’s  equations  in  component  form  become 

dy 


dEv 

dz 


=  -ico/JoJuHx 


a E,  5E 


dz  dx 


=  icojuoJuH 


dEv 

dx 


=  ia)ju0/uE[, 


and 


dHv 

dz 


=  uos0£ 


(• X)E X 


(3a) 

(3b) 

(3c) 


(4a) 


dH  dH.  .  (  w 

dz  dx  0  v  ’  y 


(4b) 


8Hv 

_ y_ 

dx 


=  -i(O£0£^x)Ez 


(4c) 


oZ _ /L 

FIGURE  1 .  Parallel  Conducting  Plates  at  x  =  0  and  x  =  a. 
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Taking  the  derivative  with  respect  to  z  of  Equation  3a  and  the  derivative  with  respect 
to  x  of  Equation  3c  and  substituting  the  result  into  Equation  4b,  Equation  4b  becomes 


d2Ey 

dx2 


k2jUs(x)-k2~\Ey 


=  0 


(5) 


It  is  assumed  throughout  that  the  total  field  components  are  given  by 

E'j'  (x,  z )  =  Ej  (x)  el{K^lut)  ;j  =  x,y,z,  and  H'°'  (x,  z)  =  (x)  el[K^w,]  ;j  =  x,y,z. 


Thus  in  the  above  components,  only  E .  (x)  and  IE  (x)  must  be  detennined.  This  in 
turn  leads  to  ordinary  differential  equations  (ODEs)  that  need  to  be  solved  as  opposed  to 
partial  differential  equations  (PDEs). 

Equation  5  is  formally  equivalent  to  a  one-dimensional  (ID)  Schrodinger  equation 
with  k_  playing  the  role  of  an  eigenvalue  and  s[x) ,  the  potential.  Once  we  have  Ey , 
from  Equations  3  a  and  3  c  we  have  II  and  Hx .  These  are  transverse  electric  (TE)  modes 
with  respect  to  z.  In  general,  for  inhomogeneous  media  it  is  not  possible  to  split  modes 
into  TE  and  transverse  magnetic  (TM)  classes,  but  when  the  permittivity  function  is 
transversely  dependent  on  a  single  variable,  it  is  possible  to  make  this  split  (References  2 
through  5). 


Returning  to  Equations  3  and  4,  we  can  obtain  an  ODE  for  Hv  also.  Taking  the 
derivative  of  Equation  4a  with  respect  to  z  and  the  derivative  of  Equation  4c  with  respect 
to  x  and  substituting  the  result  into  Equation  3b  gives 


d2H 

_ y 

dx2 


s'’(x)  dE[y 
<?(x)  dx 


\klju£{x)-k2z~]Hy=  0 


(6) 


where  Afx) 


ds^x) 

dx 


Because  we  chose  to  have  the  pennittivity  be  a  function  of  x  but  not  the 
permeability,  our  equations  for  E  and  H  are  not  symmetric.  However  from  Equations 
3  and  4,  if  Equation  6  can  be  solved  for  H  ,  then  both  Ex  and  E,  follow.  Thus  these  are 
the  TM  modes  with  respect  to  z,  the  propagation  direction. 
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3.0  EXPONENTIAL  PERMITTIVITY 


Let  the  permittivity  be  given  by 
z(x)  =bedx 


(V) 


where  b  is  often  called  the  strength  or  intensity  of  the  function  and  8  is  called  the  range. 


3.1  TECASE 

Substituting  Equation  7  into  Equation  5,  we  get 
d2E„ 


dx~ 


-  + 


~klfube^Sx  -  kl  ]  Ey  =  0 


(8) 


Solving  Equation  8,  we  obtain  for  the  general  solution  (Reference  6) 


Ey{x)  =  c  >r 


1  + 


2k 


d  2k, 
J  8 


2koy[jub  _ 


8xj  2 


+ c2r 


i- 


2k. 


d -2k. 


2K  sfjub 


Sx/2 


(9) 


The  J’s  are  Bessel  functions  of  the  first  kind.  T  denotes  the  Euler  gamma  function. 

We  are  interested  in  finite  solutions  for  a  conducting  parallel  plate  geometry  as  in 
Figure  1 .  For  now,  let  us  assume  that  one  plate  is  at  x  =  0  while  the  other  is  at  x  =  oo.  By 
doing  this,  we  see  that  the  first  solution  of  Equation  9  will  go  to  zero  as  x  — >■  oo  while  the 
second  solution  will  have  a  singularity  there. 

Thus  we  throw  away  the  second  solution  by  setting  c2=  0  in  Equation  9.  The 
electric  field  will  be  taken  as 


E y  C\  d 2kz 


2  koyfJib 


—fix/2 


8 


(10) 


co 


with  k0=  —  ,  the  free  space  wave  number,  and 


f 


Ci  =c,r 


2k 


1  + 

v  S  j 


(ID 
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Note  that  the  gamma  function  in  Equation  1 1  is  a  function  of  A:. ,  the  eigenvalue,  and 

it  will  change  for  different  modes.  To  give  Equation  10  a  dimensionless  fonn,  we 
k 

introduce  /?  =  —  where  fi  is  the  normalized  longitudinal  propagation  constant,  and  by 

K 

substitution  into  Equation  10  we  have 


Ev=c,J 


2kj^ 

5 


2  k0yfjib 


-Sx/2 


8 


(12) 


and 


Ci  =Clr 


i+ 


2 Ki 

8 


(13) 


/?  is  dimensionless,  and  thus  as  long  as  8  has  the  same  units  as  k0  (e.g.,  k{)  and  8  must 
both  have  units  of  m  1  if x  is  in  meters),  Equation  12  is  now  dimensionless. 


Now  that  E  is  known,  we  can  obtain  the  magnetic  field  components  from 
Equations  3  a  and  c 


H 


X 


(14a) 


and 


H 


z 


1  dEy 

icoju0ju  dx 


(14b) 


Thus  Hx  is  proportional  to  E 
written  as 


and  by  setting  v  = 


2  k0/3  2k0JJJb 

-J-,7  =  — ^ - > H-_  canbe 


H.  = 


Jy-i  [ve~Sx^  )  -  Jv+x  [r]e~s^2 ) 


(15) 


Equations  12,  14a,  and  15  give  the  TE  (to  z)  modes  of  this  inhomogeneous 
geometry. 
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At  this  point,  the  boundary  conditions  must  be  considered.  The  geometry  (so  far) 
consists  of  two  perfectly  conducting  parallel  plates  at  x  =  0  and  x  =  oo.  The  boundary 
conditions  are  tangential  E  is  zero  on  x  =  0,  oo  and  normal  B  is  zero  on  x  =  0,  oo.  Also 
dH 

- -  =  0  on  x  =  0,  oo  where  n  refers  to  the  nonnal  to  the  plates.  These  imply 

dn 


dH 

Ev  =  0,  H  =  0,  and - -  =  0  on  x  =  0,  oo 

dx 


(16) 


Since  E  and  Hx  are  proportional,  if  one  of  them  satisfies  its  boundary  condition, 
the  other  will  also.  Bear  in  mind  that 


dH.  _  1  d2Ey 

dx  icOjUoJu  dx1 


(17) 


We  see  from  Equation  8,  if  E  =  0  at  some  point,  - -f-  =  0  there  also.  Thus  all 

dx 

three  boundary  conditions  in  Equation  16  are  satisfied  when  E  =  0  is  satisfied. 


Now  E  must  be  0  at  x  =  0  and  at  x  =  oo.  Using  Equation  12,  E  is  zero  at  infinity 
automatically  without  restrictions  on  either  the  argument  or  the  order  parameters.  Thus 
from  Equation  12,  we  are  left  with  the  boundary  condition  (BC)  at  x  =  0: 


^1  —  J  2k :0/3 


2k0Jub 


5 


=  0 


(18) 


where  k0  = 


and  /0  is  the  linear  frequency. 


For  now  it  is  assumed  that  k0,  8  ,  b,  and  //  are  all  real  and  positive.  This  restriction 
can  be  relaxed  later.  In  this  case,  (3  is  not  only  the  propagation  constant,  it  is  also  the 
eigenvalue  of  Equation  18.  So  the  value(s)  of  [3  that  causes  a  zero  of  the  Bessel  function 
in  Equation  18  is  found  given  the  frequency  as  well  as  b,  // ,  and  8  . 

In  order  for  Equation  1 8  to  have  a  real  solution,  f3 ,  at  all,  it  is  necessary  that  the 
strength  of  the  dielectric  function  be  greater  than  a  certain  limit.  This  limit  results  from 
the  fact  that  the  first  zero  of  Jv  (if)  moves  to  smaller  values  of  77  when  v  is  decreased 
(Reference  7).  When  /?  =  0 ,  the  order  of  the  Bessel  function  in  Equation  18  is  also  zero, 
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and  the  lowest  zero  of  the  zeroth  order  Bessel  function  is  2.4048.  Thus  Equation  18  has  a 
real  solution  only  if  the  argument  of  Equation  1 8 


2  k0y[Jjb 


8 


>2.4048 


(19) 


or 


b> 


1.4457 


v^2y 


(20) 


A 


From  Equation  20,  at  lower  frequencies  b  must  be  large  and  b  can  decrease  as  the 
frequency  increases. 


3.2  TM  CASE 

For  the  TM  modes,  we  must  solve  a  slightly  more  complicated  equation  than  the  one 
used  for  the  TE  case.  Now  we  must  solve  Equation  6  when  Equation  7  is  substituted  in. 
Note  that  Equation  6  is  no  longer  in  the  fonn  of  a  ID  Schrodinger  equation  by  virtue  of 
the  first-order  derivative  in  x.  The  equation  to  solve  is 


d2H 

_ y 

dx2 


+  8- 


dH 

dx 


-  + 


~k\fibeSx  -k]\Hy=  0 


(21) 


Using  r)  =  7+N+7  as  before  and 

8 


a  = , 


2 hi 

V  s  J 


-1=77+1 


(22) 


the  general  solution  to  Equation  2 1  is  (Reference  6) 


H,  =  r-e-“p  { d[j\V-‘"2 )  +  4J_  („e^ ) 


(23) 


with 


d[  =  d1  r  (1  +  cr) ,  d^  =  d2T  (l  -  cr) 
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where  <7,  and  d2  are  general  constants.  As  previously,  we  set  d2  =  0  and  eliminate  the 
Bessel  function  of  negative  order. 

Thus  the  solution  we  will  use  is 


(24) 


At  this  point,  the  TM  case  must  proceed  a  little  differently  from  the  TE  case. 
Although  we  know  that 


E  = 


C0£Q£ 


(x) 


H„ 


(25a) 


and 


E  = 


dH„ 


-io)£0£  (x)  dx 


(25b) 


we  see  that  we  would  be  dividing  both  Hy  and  dE[y/dx  components  by  s(x)  which 
goes  to  zero  (or  at  least  becomes  very  small  at  some  point).  Thus  using  the  electric 
displacement,  D ,  at  this  point,  i.e., 


D  = 


CO£a 


-H.. 


(26a) 


and 


1  dH 
A= - z- 

-ia>£0  dx 


(26b) 


is  a  better  choice. 

Thus  Dx  is  proportional  to  H  ,  but  neither  of  these  components  has  a  boundary 
condition  associated  with  it  on  a  perfect  conductor.  D_  can  be  written  as 


D_  = 


d'Sr/ 
8  ico£n 


r/e 


-5x 


[de5*11 )  -  (rje^2 )]  +  2 e^Ja  ) 


(27) 
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Now  the  only  boundary  condition  we  must  meet  for  the  TM  modes  is  E_=  0  on  x  = 
0,  oo.  If  E_  =  0  on  these  boundaries,  then  D.=  0  on  the  boundaries  also. 

Thus  Equation  27  reduces  to  (at  jc  =  0) 

BC2=r/[ja_l(ri)-J(7+1(ri)\  +  2Jcr(ri)  =  0  (28) 

Since  cr  contains  the  eigenvalue,  J3 ,  we  solve  BC2  for  the  value(s)  of  (i  that  cause 
Equation  28  to  be  zero.  This  is  a  slightly  more  difficult  equation  to  solve  than  the  one  for 
the  TE  case. 


4.0  NUMERICAL  RESULTS 


4.1  TECASE 

For  the  TE  case,  Ey  is  the  only  nonzero  component  of  the  electric  field,  and  Hx  and 
//  are  the  corresponding  nonzero  magnetic  field  components.  Since  it  has  been  shown 
in  Section  3.0  that  all  boundary  conditions  are  satisfied,  all  numerical  results  will  focus 
on  E  alone. 

Recall  that  //  is  the  relative  (constant)  penneability,  b  is  the  relative  strength  of  the 
dielectric  function,  8  is  the  range  of  the  dielectric  function,  and  f0  is  the  linear 
frequency.  Because  8  and  k0  have  the  same  units  (i.e.,  m  1  or  cm  1 ,  etc.)  and  //  and  b 
are  dimensionless,  we  solve  Equation  18  for  /?,  the  normalized  propagation  constant. 

In  Table  1,  the  ordered  roots,  /? ,  of  Equation  18  are  given  for  8  =  1 ,  -JJJb  =  1 ,  as  a 
function  of  frequency.  /?,  is  the  eigenvalue  for  the  lowest-order  mode,  /J2  is  the 
eigenvalue  for  the  next  lowest-order  mode,  and  so  on.  As  the  frequency  increases,  more 
and  more  higher-order  modes  can  propagate  assuming  a  fixed  8  and  J~jub .  Figure  2 
shows  a  plot  of  frequency  versus  /?  also. 
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TABLE  1.  Ordered  TE  Roots,  /?,  of  Equation  18  With  8  =  1 ,  s/Jib  =  1 . 


/o  (GHz) 

A 

A 

A 

A 

A 

A 

A 

6 

0.0282 

7 

0.1206 

8 

0.1937 

9 

0.2531 

10 

0.3026 

11 

0.3446 

12 

0.3808 

13 

0.4124 

14 

0.4400 

0.0380 

15 

0.4649 

0.0796 

16 

0.4870 

0.1169 

17 

0.5070 

0.1814 

18 

0.5251 

0.1814 

19 

0.5417 

0.2088 

20 

0.5569 

0.2353 

24 

0.6069 

0.3205 

0.0914 

36 

0.6991 

0.4782 

0.3004 

0.1457 

0.0063 

94 

0.8405 

0.7223 

0.6263 

0.5421 

0.4656 

0.3948 

0.3283 
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P 

0.6  r 


0.5 


0.4 


0.3 


0.2 


0.1 


J - i - l - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - l - 1 - ' - 1 - ' - ! - 1 - 1 - 1 - 1 - * - I - 1 - 1 - 1 - 1 —  Freq(GHz) 

6  8  10  12  14  16  18  20 

FIGURE  2.  TE  Longitudinal  Propagation  Constant,  /? , 

V ersus  Frequency  (//  =  b  =  8  =  1 ) . 


Figure  2  gives  the  dispersion  curves  of  this  parallel  plate  geometry  with  an 
inhomogeneous  permittivity  of  exponential  form.  For  a  given  8 ,  ju ,  and  b,  once  the 
roots  of  Equation  1 8  are  found  at  a  given  frequency,  one  can  plot  the  normalized  electric 
field  versus  x  (cm)  (x  is  the  distance  along  the  transverse  direction  between  the  plates). 
Notice  that  so  far  we  have  still  assumed  that  the  second  plate  is  at  x  =  oo.  Figures  3(a), 

(b),  (c),  and  (d)  show  plots  of  normalized  E  versus  x  for  8  =  1 ,  yfjub  =  1 ,  and 
f0  =  24  GHz .  At  24  GHz,  there  are  three  real  roots  of  Equation  18  (see  Table  1).  Using 
/?  =  /?,  in  Ey  gives  Figure  3(a),  the  lowest-order  mode;  using  J3  =  /?,  in  E  gives 
Figure  3(b),  the  second  mode;  and  /?  =  /?3  in  E  gives  Figures  3(c)  and  3(d),  the  third 

mode.  Figure  3(d)  is  the  same  plot  as  Figure  3(c),  only  on  an  expanded  scale.  In 
Figure  3(a),  note  that  E  is  highly  concentrated  near  the  conducting  plate  at  x  =  0  and 

trails  off  to  zero  as  x  increases.  (The  asymptotic  forms  of  the  fields  are  given  in 
Appendix  A.)  The  orthogonality  property  of  different  propagating  modes  at  a  given 
frequency  is  shown  in  Appendix  B.  Also  the  Poynting  vector  and  the  power  flow  along 
the  inhomogeneous  waveguide  are  derived  in  Appendix  C. 
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Ey  (normalized) 


x  (cm) 


FIGURE  3.  Normalized  Electric  Field,  Ey  ,  Versus  x 
for  the  TE  Case  (//  =  b  =  8  =  1,  /0  =  24  GHz) . 
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Ey  (normalized) 


(d)  Third  lowest-order  mode  on  larger  scale. 
FIGURE  3.  (Contd.) 
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Ev  is  not  exactly  zero  unless  x  =  oo,  but  we  do  not  need  to  require  this.  We  see  that 

under  the  proper  conditions,  there  is  actually  no  need  for  the  second  conducting  plate  at 
all.  The  second  boundary  of  the  inhomogeneous  region  can  be  set  at  a  value  of  x  such 
that  the  magnitude  of  the  field  E  is  less  than  some  prescribed  amount,  say  E  <  1(T6  at 

x  =  a’ . 


Thus  using  Equations  10  and  18,  after  5 ,  b,  //  and  f0  are  assumed,  one  substitutes 
into  Equation  10  (with  the  appropriate  (i  from  BC\  )  giving 


(2k  I — 
Ey(x  =  a')  =  J 2_w  -tV/Ai 
5  V 


—da' 


5 


(29) 


By  enforcing  a  prescribed  value  such  as  Ey  <  10  6 ,  we  can  calculate  a  value  of  a' 
(this  is  not  the  position  of  the  second  conducting  plate).  The  second  plate  has  been 
removed,  and  a'  is  the  distance  at  which  the  field  is  appropriately  small. 

Figure  4  shows  a  plot  of  the  nonnalized  electric  field  (lowest-order  mode)  versus  jc 
for  a  higher  frequency,  f0  =  35  GHz  and  ju  =  8  =  b  =  1 .  As  frequency  increases,  the  form 
of  E  becomes  narrower  and  more  closely  confined  to  the  plate  at  x  =  0.  As  frequency 
increases,  E  <  10  6  will  be  enforced  by  using  smaller  and  smaller  values  of  a  . 
However,  as  the  higher-order  modes  are  considered  (as  in  Figure  3),  the  higher  the  mode 
at  a  given  frequency,  the  less  tightly  it  is  confined  to  jc  =  0  (see  Figures  3(c)  and  3(d)). 


Figure  5  indicates  what  happens  when  ju  =  \  =  S  but  b  is  increased  to  4  at 
f0=  35  GHz  .  By  increasing  either  //  or  b,  one  can  confine  the  field  even  more  tightly  to 
x  =  0  (note  that  Figure  5  only  runs  to  jc  =  2.0,  whereas  Figures  3  and  4  run  to  jc  =  4.0). 
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Ey  (normalized) 


FIGURE  4.  Normalized  Electric  Field,  Ey  ,  Versus  x  for  the  TE  Case 
[ju  =  b  =  5  =  1,  /0  =  35  GHz)  (Lowest-Order  Mode). 


Ey  (normalized) 


FIGURE  5.  Normalized  Electric  Field,  Ey  ,  Versus  x  for  the  TE  Case 
[ju  =  S  =  \,  b  =  4,  /0  =  35  GHz)  (Lowest-Order  Mode). 
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The  cutoff  frequency  for  each  mode  occurs  at  (5  =  0  or  from  Equation  19, 


Jojc8 

471-J-jub 


l  =  1,2,3... 


(30) 


where  j0 1  are  the  zeros  of  the  zeroth  order  Bessel  function.  Thus  (see  Figure  2),  the 
cutoff  frequency  of  the  lowest-order  mode  is  given  by  the  lowest  zero  of  J()  ( j0  l )  (and 
with  8  =  1  =  ju  =  b  ) 


fjo  =5.74  GHz 
/W  =13.18  GHz 

and  so  on  for  the  higher-order  modes. 

4.2  TM  CASE 

For  the  TM  case,  H  is  the  only  nonzero  component  of  the  magnetic  field,  and  E 
and  E.  are  the  corresponding  nonzero  electric  field  components.  However,  in  this  case 
there  is  no  boundary  condition  that  Hv  (a  tangential  component)  must  obey.  Thus,  we 
work  with  E z  since  E,  must  be  zero  on  both  parallel  plates.  (As  discussed  in 
Section  3.0,  we  use  Dz).  In  this  case,  Equation  28  is  solved  for  /?,  the  nonnalized 
propagation  constant  for  the  TM  modes.  In  Table  2,  the  ordered  TM  roots,  /? ,  of 
Equation  28  are  given  for  <7=1,  J~jub  =  1 ,  as  a  function  of  frequency.  As  the  frequency 
increases,  more  and  more  higher-order  modes  can  propagate  assuming  a  fixed  8  and 
JJib  .  Figure  6  is  a  plot  of  frequency  versus  J3  for  the  TM  case.  In  general,  the  [i ’s  for 
the  TM  case  are  higher  than  those  for  the  TE  case  at  the  same  frequency.  Figures  7(a), 
(b),  and  (c)  show  a  plot  of  D_  (normalized)  versus  x  in  cm  for  the  first  three  TM  modes 

with  8  =  ju  =  b  =  1  and  fQ  =  24  GHz .  The  normalized  D_  component  is  also  tightly 

confined  to  the  conducting  plate  at  x  =  0  and  goes  to  0  as  x  increases,  again  holding  the 
TM  fields  within  the  inhomogeneous  region. 
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TABLE  2.  Ordered  TM  Roots,  /? ,  of  Equation  28  With  8  =  1 ,  ^J~Jub  =  1 . 


/o  (GHz) 

A 

A 

A 

A 

A 

A 

A 

6 

0.1788 

7 

0.3678 

8 

0.4639 

9 

0.5283 

10 

0.5758 

11 

0.6128 

12 

0.6426 

13 

0.6674 

14 

0.6882 

0.1281 

15 

0.7061 

0.1907 

16 

0.7217 

0.2367 

17 

0.7354 

0.2743 

18 

0.7475 

0.3064 

19 

0.7584 

0.3345 

20 

0.7682 

0.3595 

24 

0.7993 

0.4384 

0.1706 

36 

0.8526 

0.5756 

0.3782 

0.2089 

0.0302 

94 

0.9263 

0.7773 

0.6719 

0.5823 

0.5022 

0.4286 

0.3599 
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0.6 


0.4 


0.2 


Freq  (GHz) 


FIGURE  6.  TM  Longitudinal  Propagation  Constant,  (3 , 
Versus  Frequency  (//  =  b  =  8  =  1 ) . 


Dz  (Normalized) 


FIGURE  7.  Normalized  Electric  Displacement,  D_  Versus  x. 


Figures  7(a),  (b),  and  (c)  are  the  three  lowest-order  TM  modes 
^  =  b  =  S  =  \,f0  =  24  GHz). 
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Dz  (Normalized) 


Dz  (Normalized) 


FIGURE  7.  (Contd.) 
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The  TM  mode  cutoff  frequencies  are  found  when  (3  =  0  or  cr  =  1 .  Using  a  =  1  in 
Equation  28,  Equation  28  becomes 

?7[‘/o(?7)-‘/2(77)]  +  2*/i(77)  =  0  (31) 

,  .  2J,(tj)  .  . 

Writing  J2yrj)  = - - — -  ,  Equation  3 1  gives 


JM  =  0  (32) 

as  the  cutoff  condition.  Thus  the  TM  mode  cutoff  condition  is  the  same  as  that  for  the  TE 
modes  given  by  Equation  30.  (This  can  also  be  seen  numerically  from  Figures  2  and  6.) 


5.0  COMPARISON  OF  INHOMOGENEOUS  AND  HOMOGENEOUS 
PARALLEL  PLATE  WAVEGUIDES 


5.1  TECASE 

It  is  very  interesting  to  compare  the  simple  parallel  plate  waveguide  containing  a 
homogeneous  medium  and  the  parallel  plate  waveguide  containing  a  medium  with 
inhomogeneous  permittivity  in  exponential  form.  The  TE  electric  field  in  the 
homogeneous  case  is  given  by  Reference  1 


Ev  =  E0  sin  eiKze~M ;  m  =  1,2,3... 


(33) 


for  the  TE  modes,  while  in  the  inhomogeneous  case,  it  is  given  by  Equation  10.  E 
given  by  Equation  33  is  a  sinusoid  that  is  periodic.  E  is  zero  at  x  =  0,  but  there  must  be 
a  second  conductor  at  x  =  a  in  order  to  have  the  field  go  to  zero  there.  The  higher  order 
modes  split  the  interval  fromx  =  0  to  a  into  equal  subintervals.  An  infinite  series  of  these 
guides  could  be  placed  along  x  =  0,  a\  x  =  a  to  2a;  x  =  2a  to  3 a,  and  so  on.  The 
dispersion  relation  for  the  homogeneous  parallel  plate  waveguide  is  given  by 


k:=klsju- 


/  \2 

'  ni7t  ' 


V  a  J 


(34) 
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A  plot  of  k_  versus  k0  is  shown  in  Reference  1  on  page  135.  It  is  pointed  out  that 


for  the  mth  mode  at  some  point,  k  will  become  imaginary  if  k0  < 


1)171 


.  The  wave 


then  becomes  evanescent  and  attenuates  exponentially  in  the  longitudinal  direction. 
Compare  this  with  the  dispersion  relation  for  the  inhomogeneous  case 


J 


2 M 
s  V 


Ik, 

5 


]=o 


(35) 


{K=KP) 


The  homogeneous  dispersion  relation  is  explicitly  known  while  the  one  above 
(Equation  35)  is  implicit.  But  while  k,  in  Equation  34  can  become  purely  imaginary  for 

J717T 

k0  <  — j= ,  if  /?  is  assumed  to  be  imaginary  in  Equation  35,  consider  what  happens. 
a-sjs/j 


Let  (3  =  i% ,  x  e  D  and  then  the  dispersion  relation  becomes 


J 


2*o  <Z 


6 


=  0 


(36) 


But  this  equation  is  complex,  and  while  its  real  and  imaginary  parts  do  have  zeros, 
they  are  not  the  same  zeros.  Thus  there  are  no  purely  imaginary  solutions  ix  to  the 
above  (for  ju ,  8  ,b  positive  and  real).  For  the  inhomogeneous  parallel  plate  waveguide, 
there  are  no  evanescent  waves  that  purely  attenuate  along  z  when  // ,  8 ,  and  b  are 
assumed  to  be  positive  and  real.  When  both  the  order  and  the  argument  of  Equation  35 
are  assumed  to  be  complex,  complex  roots  will  occur  (References  8  and  9).  This  is 
interesting  in  that  for  the  TE  case,  the  homogeneous  guide  containing  a  lossless  isotropic 
medium  has  a  denumerably  infinite  set  of  modes  at  any  given  frequency,  most  of  which 
are  evanescent  and  thus  cut  off  with  a  few  (or  perhaps  only  one)  propagating  modes.  The 
inhomogeneous  waveguide  containing  a  lossless  isotropic  inhomogeneous  medium  has  a 
spectrum  at  any  given  frequency  that  is  completely  real  and  finite.  At  a  given  frequency 
(above  the  first  cutoff)  one  has  a  finite  number  of  propagating  modes.  Purely  evanescent 
modes  do  not  exist. 
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For  the  homogeneous  parallel  plate  waveguide,  the  cutoff  frequency  is  (References  1 
and  10) 


me 


2a^[sju 


m  =  1,2,3,... 


(37) 


where  c  is  the  speed  of  light.  There  are  two  ways  to  lower  this  cutoff  frequency.  One  is 
by  increasing  e  and  // ,  and  the  other  is  by  making  a ,  the  plate  separation,  larger. 
Similarly,  when  considering  Equation  30,  which  gives  the  cutoff  frequencies  of  the 
inhomogeneous  parallel  plate  waveguide,  the  only  way  to  lower  the  inhomogeneous 
cutoff  frequencies  is  to  either  increase  //  and  b  (analogous  to  increasing  //  and  s  in 
Equation  37)  or  to  decrease  the  value  of  8  which  governs  the  range  of  the  permittivity. 
From  Equation  29,  decreasing  8  implies  that  a'  must  be  increased  to  give  the  prescribed 
value  of  E  <  1CT6  at  x  =  a' . 

y 


5.2  TM  CASE 

In  the  TM  case,  the  homogeneous  parallel  plate  waveguide  has  a  D_  component 
given  by  Reference  1 

D  =Z^Asin  kxxei{Kz-M)  (38) 

co 


TY17Z 

where  kx  = -  (just  as  in  the  TE  case  except  that  m  can  be  zero  in  the  TM  case).  The 

a 

D  component  in  the  inhomogeneous  case  is  given  by  Equation  27.  The  TM  dispersion 

relation  for  the  homogeneous  parallel  plate  waveguide  is  the  same  as  for  the  TE  case. 
The  TM  dispersion  relation  for  the  inhomogeneous  guide  is  given  by  Equation  28. 


Again,  as  with  the  TE  case,  k,  in  Equation  34  can  become  purely  imaginary 


tor 


kn  < 


nm 


.  We  look  at  Equation  28  to  see  what  would  happen  if  /?  were  assumed  to  be 


a^fsju 

purely  imaginary.  For  /?  =  i% ,  %  e  □ 


a  = 


iX 

8  j 


+  1 


(39) 


or 
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a  = 


l  S  J 


2 


(40) 


2k  y 

Thus,  while  a  will  still  be  real  if  — —  <  1 ,  substitution  of  /3  =  i%  into  Equation  28 

8 

does  not  result  in  real  solutions  for  % .  Although  we  have  only  considered  a  limited 

2k  y 

number  of  cases,  no  roots,  % ,  have  been  found  that  satisfy  Equation  28.  For  — - —  >  1 , 

8 

<j  becomes  imaginary  and  again  no  roots,  % ,  have  been  found  that  satisfy  Equation  28 
for  /j ,  b,  8  real  and  positive.  (Thus,  the  inhomogeneous  parallel  plate  waveguide 
spectrum  again  is  completely  real  and  finite,  and  there  are  no  purely  evanescent  TM 
modes.)  If  any  of  ju ,  b,  or  8  is  allowed  to  be  complex,  again  complex  roots  can  occur. 


6.0  CONCLUSION 


A  parallel  plate  waveguide  containing  a  region  with  a  transverse  exponential 
permittivity  function  of  position  has  been  analyzed.  It  has  been  shown  that  it  is  possible 
to  eliminate  one  of  the  waveguide  conducting  plates  and  still  keep  the  electromagnetic 
field  closely  confined  to  the  remaining  plate  and  held  within  the  inhomogeneous  region. 
A  comparison  between  the  inhomogeneous  parallel  plate  geometry  and  the  well-known 
homogeneous  parallel  plate  waveguide  shows  a  number  of  interesting  differences. 

First,  for  lossless  media  (i.e.,  ju  ,  8 ,  b  positive  and  real),  the  eigenvalue  spectrum  of 
the  inhomogeneous  parallel  plate  waveguide  is  real  and  finite  at  a  given  frequency  as 
opposed  to  the  spectrum  of  the  homogeneous  waveguide,  which  is  denumerably  infinite 
at  a  given  frequency  and  composed  of  eigenvalues  that  are  either  purely  real  (propagating 
modes)  or  purely  imaginary  (attenuated  modes).  Thus,  the  inhomogeneous  case  has  no 
purely  attenuated  modes. 

Second,  the  eigenfunctions  of  the  inhomogeneous  case  are  not  periodic  sinusoids  but 
are  Bessel  functions  containing  an  exponential  argument  which  decrease  in  amplitude  as 
the  transverse  distance  increases.  While  the  eigenfunctions  of  the  homogeneous  parallel 
plate  waveguide  divide  the  interval  between  the  plates  into  equidistant  sections,  the 
eigenfunctions  of  the  inhomogeneous  case  can  have  nodes  at  any  point  along  the  interval 
from  x  =  0  to  x  =  a  ,  and  the  number  of  nodes  is  small.  All  fields  in  the  transverse 
direction  attenuate  exponentially.  This  sharp  amplitude  decrease  allows  the  removal  of 
the  second  conducting  plate  and  causes  the  field  components  to  be  closely  confined  to  the 
remaining  plate. 
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Third,  the  cutoff  frequencies  for  the  inhomogeneous  case  are  dependent  on  the  zeros 
of  the  zeroth  order  Bessel  function,  while  the  cutoff  frequencies  for  the  homogeneous 
case  are  dependent  on  the  set  of  positive  integers. 

All  numerical  results  have  been  based  on  the  assumption  of  lossless  media.  In  future 
work,  lossy  media  (//  and  b  complex)  will  be  considered.  Also,  8  complex  will  be 
assumed,  which  changes  the  purely  exponentially  attenuating  pennittivity  function  into 
one  that  is  sinusoidally  oscillating  and  damped. 
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Appendix  A 

ASYMPTOTIC  FORM  OF  Ey  AND  //, 


A-l 
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For  the  TE  case,  E  is  given  by  Equations  10  and  11.  For  large  x,  the  argument  in 
Equation  10  becomes  very  small,  and  Jv  [tje  ‘ix/2  j  can  be  approximated  by  the  first  term 
of  its  series  expansion.  Thus  for  large  x, 


Ey(x)*C 1 


rTjeSxl2^ 


or 


Ey(X) 


r  7j^ 

v2. 


-Svxj  2 


The  total  field  becomes 


p  total 


'M' 

v  s  J 


Jr k0/3x+ik0/3z-io>t ) 


(A-l) 


(A-2) 


(A-3) 


where  (3 ,  77 ,  v  are  assumed  to  be  both  positive  and  real. 


For  the  TM  case,  for  large  x, 


Hy  (x)  *  dx 


kJ3_ 

8  j 


(A-4) 


where  a  =  \lv2  +  1  and  i] ,  and  cr  are  given  by  Equation  22  and  previously.  The  total 
asymptotic  field 


tt  total 

Hy 


(A-5) 


For  large  x,  all  field  components  fall  off  as  exponentials. 
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Appendix  B 

ORTHOGONALITY  OF  FIELD  COMPONENTS  FOR 
THE  INHOMOGENEOUS  PERMITTIVITY  4.v)  =  be& 
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In  the  TE  case,  the  electric  field  has  only  a  y  component, 

Ey=c;jv(r,e-Sxl2)  (B-l) 

with  v  and  77  given  in  Section  3.0. 

For  large  enough  frequencies,  it  is  possible  to  have  more  than  one  p  value 
satisfying  the  boundary  condition  at  x  =  0  given  by  Equation  18.  When  k0 ,  b,  // ,  8  are 
real,  positive,  and  have  fixed  values,  if  there  are  two  solutions  to  Equation  18  (i.e.,  /?, 
and  P2 ),  then  the  corresponding  eigenfunctions  are 


with  v1 


2k„Jbfi 

s  ,n  s 


and 


(B-2a) 


with  v2  = 

2  8 

E (1*  and 

y  y 


77  as  above, 
are  orthogonal  if 


00 

J E J1*  (x)w(x)dx  =  0 

0 

where  w(x)  is  some  weighting  function. 
For  now  we  set  w(x)  =  1 . 

Let 


(B-2b) 


(B-3) 


I  =  \ (x)  Ef  (x) dx  =  c|,(lV(2,jX  (r/e~dx/1 )  JVi  (r/e~Sx/2 )  dx 


(B-4) 
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Making  the  variable  change, 

-Sxfl 


u~e 


(B-5) 


we  obtain 


dx  =  - 


2  du 
8  u 


(B-6) 


Substitution  of  Equations  B-5  and  B-6  into  Equation  B-4  gives 


1  = 


r(i)c  '(2)  r7,,  (nu)Jt,2  (nu)du 
1  J  n 


(B-7) 


Upon  integration,  /  becomes 


/  = 


'  'W  '6) 
.  C1  C1 


(m/7)U2  {uri)~Jv,  (ut))[ut)JV2_,  (urj)  +  (v,  -v2)JVi  (ur/)_ 


(B-8) 


Since  J  (7)  =  0  ,  and  Jv  (//)  =  0  from  the  boundary  condition,  upon  evaluation  at 
the  limits,  /  reduces  to  zero.  Thus,  two  eigenfunctions  with  different  ft  ’s  at  the  same 
frequency  are  orthogonal.  The  TM  case  is  more  complicated  but  done  similarly. 

If  the  orders  of  the  eigenfunctions  are  different  and  the  arguments  are  different  also, 
orthogonality  does  not  hold. 
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Appendix  C 

THE  POYNTING  VECTOR  AND  POWER  FLOW 


C-l 


NAWCWD  TP  8738 


This  page  intentionally  left  blank. 


C-2 


NAWCWD  TP  8738 


TE  CASE 


The  complex  Poynting  vector  is 

5  =ExH* 

(C-l) 

=  Eyyx[H*xx  +  H*j] 

(C-2) 

=  -EyH*xz  +  EyH*zx 

(C-3) 

where  *  implies  complex  conjugate.  Ey  is  given  by  Equation  12,  Hx  and  H.  by 
Equation  14.  With  n  ,  b,  8 ,  k0,  and  (5  all  real  and  positive 

E,=c;jy(, (C-4) 

«;=—  c[j\nes^)  (C-5) 

and  H*z  is  given  by  the  complex  conjugate  of  Equation  15.  (The  z  and  time  dependencies 
have  cancelled.) 

From  Equations  C-3,  C-4,  and  C-5,  we  see  that  the  z  component  of  the  Poynting 
vector  is  real,  i.e., 

S.  =  ^_KP_j2  i-sxn  \  (C-6) 

O)jU0jU 

while  the  x  component  of  the  Poynting  vector  is  imaginary. 

Because  the  power  is  defined  as 

P  =  Re  JJ  (ExH*)ndz  (C-7) 

cross  section 

the  power  flowing  in  the  z  -direction  is 

f2  7  o  00 

P  =  f  J2v  ( neSxl1 )  dx  (C-8) 

O)fU0fU  J0 
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or  using  the  same  variable  substitution  as  in  Appendix  B, 


P  = 


2c;\p\J2v{m) 


Sa>jU0jU 


f 


du 


(C-9) 


Upon  integration  the  power  becomes 


P  = 


c[%P 


Scoju0ju  v[r(l  +  v)] 


7 


n2k 


F 

2 1  3 


{n,^-  +  v},  {1  +  v,  1  +  v,  1  +  2v} ,  - rf 


(C-10) 


where  2F3  is  a  generalized  hypergeometric  function  with  two  numeratorial  parameters 
and  three  denominatorial  ones  as  defined  in  Mathematica. 


TM  CASE 


In  this  case,  the  complex  Poynting  vector  is 


S  =  ExHyz  -  EzHyx 


(C-ll) 


Assuming  all  variables  real  and  positive,  H*v  is  given  by  Equation  24  and  Ex  and 
E.  by  Equation  25.  The  z  component  of  the  Poynting  vector  is  real  and  given  by 


= 


KP 


a)£0e(x) 


V2  j 


d;2e-SxJ2a(rje-Sxl2) 


(C-12) 


or 


= 


KP 


CO£.b 


V 


(C-13) 


Using  Equation  C-7  and  the  same  procedure  as  in  the  TE  case,  the  power  for  a  TM 
mode  is 


P  = 


KK1  f  ;7v 


CQ£ob  \2J  o-^r(l  +  cr  )T 


2  |[j-j  2F2  |o-,^-  +  <t|,{1  +  cj,1  +  o-,1  +  2o-},-72 


(C-14) 
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